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Accurate parameter estimation of gravitational waves from coalescing compact binary sources is 
a key requirement for gravitational-wave astronomy. Evaluating the posterior probability density 
function of the binary's parameters (component masses, sky location, distance, etc.) requires com- 
puting millions of waveforms. The computational expense of parameter estimation is dominated 
by waveform generation and scales linearly with the waveform computational cost. Previous work 
showed that gravitational waveforms from non-spinning compact binary sources are amenable to a 
truncated singular value decomposition, which allows them to be reconstructed via interpolation at 
fixed computational cost. However, the accuracy requirement for parameter estimation is typically 
higher than for searches, so it is crucial to ascertain that interpolation does not lead to significant 
errors. Here we provide a proof of principle to show that interpolated waveforms can be used to 
recover posterior probability density functions with negligible loss in accuracy with respect to non- 
interpolated waveforms. This technique has the potential to significantly increase the efficiency of 
parameter estimation. 



I. INTRODUCTION 

Astronomy and tests of fundamental physics with grav- 
itational waves from compact binary coalescence (CBC) 
will ultimately be limited by our ability to estimate the 
binary's source parameters from the gravitational-wave 
signal [e.g.,[TH3]. CBC sources with total masses in the 
range 2 M Q < Mr < 500 M Q will be amongst the prime 
sources for advanced LIGO [I] and VIRGO [5] when they 
begin operating around 2015 [6]. 

In a Bayesian treatment of parameter estimation, one 
is interested in the probability distribution of the set of 
source parameters of the underlying model given an ob- 
served stretch of data. Waveform computation represents 
the majority of the computation cost in the Bayesian 
analysis of CBC sources, so the total computational cost 
scales roughly linearly with waveform generation. This 
becomes burdensome when one needs to explore a large 
dimensional parameter space as the number of waveform 
computations is large, e.g., O(10 6 ) [7]. 

Recently, Cannon et al. [8] showed that a truncated 
singular value decomposition (SVD) can be applied to 
gravitational-wave template banks which span the two 
parameters describing the masses of the coalescing bi- 
nary. The SVD decomposes the bank into a set of "ba- 
sis templates" and projection coefficients. In general, 
the number of basis templates is much less than the 
total number of templates in the bank. Furthermore, 
the projection coefficients can be interpolated across the 
domain of the bank [pj. Template waveforms can thus 
be interpolated. Because the intrinsic parameter space 
of CBC sources with non-spinning components is two- 
dimensional (the two mass parameters), it is possible 
to set up the waveform computation for parameter es- 



timation such that the waveform calculations are done 
by interpolation alone. However, the errors incurred 
from interpolation could, in principle, affect parameter- 
estimation accuracy. 

In this paper, we describe the application of SVD- 
intcrpolatcd waveforms to CBC parameter estimation. 
For a simulated data set containing a gravitational 
wave signal we provide a proof of principle that SVD- 
interpolated waveforms can be used for parameter es- 
timation without significantly affecting the accuracy of 
the inferred probability distributions of the source pa- 
rameters. We further show that the computational cost 
of using interpolated waveforms is around an order of 
magnitude less than that of commonly-used time-domain 
waveform families. This technique has the potential to 
increase the computational efficiency of CBC parameter 
estimation when the computational cost is dominated by 
waveform computation. Our application of the SVD is 
limited to a small patch of parameter space about the 
injected signal value. 

This paper is organized as follows. In Sections |Tl| and 
|III| we outline the principles of parameter estimation for 
CBCs and interpolating template waveforms based on 
the SVD, respectively. In Sec. |IV| we describe the appli- 
cation of SVD-interpolated waveforms to parameter esti- 
mation. In Sec. [V] we compare the results of using inter- 
polated and non-interpolated waveforms for parameter 
estimation and compare the computational cost of inter- 
polation to using non-interpolated waveform families. In 



Sec. VI we consider the future of using SVD-intcrpolatcd 
waveforms for parameter estimation and discuss the tech- 
nical requirements of implementing these waveforms in 
parameter-estimation pipelines. 
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II. CBC PARAMETER ESTIMATION WITH 
MCMC 

The central quantity of interest in Bayesian parameter 
estimation is the posterior probability density function 
(PDF) of a set of parameters 9 which parameterize a 
model, H, assumed to describe a data set d. The PDF 
is related to the likelihood function and prior probability 
via Bayes' theorem and is given by 



p(0\d,U) = 



V(9\H) £(d\0,H) 
p(d\H) 



(2.1) 



where £(d\9,H) is the likelihood function and V(0\T~L) is 
the prior probability which encodes our a priori belief 
in the distribution of 9. The quantity in the denomina- 
tor, p(d\H), is known as the evidence and is an overall 
normalization factor which we will not deal with here. 

The CBC parameter vector is high-dimensional. 
The phasing and amplitude of a waveform from a 
non-spinning coalescing compact binary source is con- 
trolled by two mass parameters, the chirp mass Ai = 
(TOim 2 ) 3 ^ 5 /(m 1 + m 2 ) 1 ^ 5 an d symmetric mass ratio r\ — 
(ra\m-i) I (m\ + m-i) 2 , where ra\ and rri2 are the com- 
ponent masses of the binary. In addition, a gravitational 
wave source with respect to the earth is specified by loca- 
tion dependent parameters. These are the distance from 
the Earth D, inclination t, right ascension a, declination 
S, polarization phase tp and time and phase at coales- 
cence, t c and <p c . In general, the CBC parameter vec- 
tor 9 is nine-dimensional for circular binaries with non- 
spinning components. 

One of the goals of gravitational- wave astronomy is to 
estimate the PDF of the parameters of a candidate gravi- 
tational wave source in order to assign a meaningful prob- 
ability to our measurements of the source properties and 
demographics. To compute the right hand side of (2.1 1, 



we directly evaluate the likelihood function, C(d\9,%). 
Under the hypothesis that the data, d, consists of Gaus- 
sian, stationary noise n and a gravitational-wave signal 
h(9), the likelihood is a Gaussian [10] : 



£(d\9,H) oc e 



-(d-h(9)\d-h{6))/2 



(2.2) 



where (a\b) is the usual noise- weighted inner-product, 



(o|6) = 45R 



df 



In 



d(f)b*(f) 

Sn(f) ' 



(2.3) 



and S n (f) is the detector's noise power spectral density 
(PSD). The limits of integration correspond to the band- 
width of the detector. A significant computational cost 
of evaluating the likelihood comes from computing the 
template waveform h{9) at each point in the parameter 
space. 

The full PDF is multi-dimensional and to get estimates 
on individual parameters we work with the marginalized 



PDF of either a single parmeter 9a G 9 or pairs of param- 
eters (9a, 9b) & 9. Writing 9 = (9a, ©), the marginalized 
one-dimensional PDF of 9a is thus: 



p(9 A \d,U) 



dQp(9\d,H), 



(2.4) 



and the extension to pairs is trivial. 

To efficiently evaluate the likelihood function we typ- 
ically use a stochastic sampling algorithm. Here we em- 
ploy Markov-chain Monte Carlo (MCMC), whose appli- 
cation to gravitational-wave parameter estimation is de- 
scribed in [TTj. The basic element of a MCMC is the 
"Markov chain" which directly samples the posterior dis- 
tribution, after an initial "burn-in" is discarded. 

We use the stationary phase approximation (SPA) in- 
spiral waveforms for both our simulated signal and tem- 
plate model. This allows us to directly "inject" the signal 
waveform into simulated frequency-domain noise without 
performing an additional Fourier transformation, which 
could introduce spurious artifacts related to the abrupt 
in-band termination of the time-domain waveform. For 
fixed a, S, l, and (p, the post-Newtonian frequency- 
domain waveform has the form 



Hf) = 



A(M, V ;f) 
D 



,2Trift c -i<p c +i1>(M,ri;f) 



(2.5) 



where expressions for the amplitude A(Ai,r);f) and 
phase ^(Ai, rj; /) can be found in for TaylorF2 post- 
Newtonian waveforms at zeroth order in amplitude and 
3.5 post-Newtonian order in phase. 

In general, the template model will be determined by 
the needs of a particular analysis, however for this proof 
of principle we will not consider time-domain template 
waveforms. 



III. INTERPOLATING TEMPLATE 
WAVEFORMS USING THE SVD 

The interpolation scheme used in [9] is based on the 
truncated SVD of a gravitational-wave template bank. 
The SVD decomposes the template bank into a set of or- 
thogonal basis templates, the number of which equals the 
number of templates in the bank, and projection coeffi- 
cients. Any template in the bank can be reconstructed 
from the bases weighted by the appropriate projection 
coefficients. However not all the bases are required to ap- 
proximately reconstruct the waveforms. Truncating the 
SVD reduces the number of unique basis templates; we 
choose to truncate so that the norm of any reconstructed 
template is conserved to a level of ~ 10 -5 (c.f. Eqs. (27) 
and (28) in [8]). The coefficients can be interpolated 
within the domain of the template bank which takes us 
from a discrete description of the template bank to a con- 
tinuous one. Any template waveform in this domain is 
then approximately recovered by a linear combination of 
the basis templates weighted by the appropriate interpo- 
lated projection coefficients. 
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Specifically, the SVD of a template bank of gravita- 
tional waveforms allows them to be written as a linear 
combination of basis waveforms with projection coef- 
ficients M^Mki^i) where the indices k and I enumerate 
a particular template in the bank. This implicitly as- 
sumes that the template bank follows a rectangular grid 
in (M , rj) . The waveform can thus be written 



h(M k ,m) = J2 M ^ M 



(3.1) 



To find the basis vectors we compute the SVD of 
a real- valued N x L matrix, with N < L, H whose rows 
are template waveforms. The SVD factors H as 



N 



(3.2) 



where v is a unitary matrix of reconstruction coeffi- 
cients, a is a vector of singular values and u is a ma- 
trix of orthonormal basis vectors u. In the above for- 
malism the projection coefficients are related to the 
reconstruction-coefficient matrix and vector of singular 

values by M^(M k ,m) : = o>( v (2^)m + *"(2i/+i)/*) where v 
enumerates a (complex) waveform in the template bank. 
The indices (k, I) are related to v by the specific packing 
of the matrix H, the details of which are irrelevant. 

The projection coefficients can be interpolated and we 
employ the method of [5], using Chebyshev polynomials 
of the first kind (c.f. Eqs. (7) and (8) in [5]). An "interpo- 
lat cd w aveform", h' , can thus be constructed according 
to (3.1 ) from a linear combination of interpolated coeffi- 



cients and basis waveforms: 



(3.3) 



This representation of h! is continuous and hence any 
waveform within the domain of the original template 
bank can be recovered by interpolation. In general the 
accuracy of interpolated waveforms depends directly on 
the density of the template bank within a particular do- 
main. Below we illustrate the application of the SVD to 
parameter estimation. 



IV. PARAMETER ESTIMATION USING 
INTERPOLATED WAVEFORMS 

We will compare the marginalized one-dimensional 
PDFs obtained by using an interpolated template wave- 
form family to those generated using the standard, non- 
interpolated template waveform family for the same data 
set. For illustration we consider a toy example with five 
free parameters. We generate a single-detector data set 
d containing a signal waveform h{6) and Gaussian sta- 
tionary noise n with a power spectral density roughly 



matching that of initial LIGO [13]. By only having five 
free parameters we effectively set the prior on the other 
four to be delta functions centered on the signal values. 
We chose to fix the sky position (a, 8) and the inclination 
and polarization angle (t, ip) of the template waveforms 
such that they arc not searched over. 

Because interpolation is carried out in the mass space 
only, we focus on exploring the effects of interpolation 
on mass parameters and parameters that are known to 
be very strongly correlated with masses (time and phase 
of coalescence and distance). If the accuracy of the re- 
covery of these parameters is unaffected by interpolation, 
we do not expect the angular parameters to be affected, 
either. However, it is important to realize that the abso- 
lute accuracy with which some parameters, particularly 
distance, are estimated is improved by fixing sky location 
and orientation parameters and lifting corresponding de- 
generacies. Thus, the measurement uncertainties inferred 
below should not be considered typical. Since we demand 
that systematic biases from using interpolated templates 
are smaller than statistical measurement uncertainties, 
the improvement in the accuracy of distance measure- 
ment means that we are being conservative in evaluating 
the quality of SVD-intcrpolatcd parameter estimation. 

The signal contained in the data set has 
source parameters (M., rj, D, t c , <p c ) = 
(7.45 M Q , 0.247, 33 Mpc, s, 2.16) and we have 
omitted the others which are not searched over. The 
signal has a signal-to-noise ratio SNR = 14.8. The 
frequency-domain data is sampled at A/ = 1/32 Hz 
with a maximum frequency of 512 Hz. 

The prior distributions are set as follows. We use 
a uniform prior on log_D and r\ with ranges D G 
[1 Mpc, 100 Mpc] and r] G [0.175, 0.250]. We use a prior 
on M of the form V{M\H) oc M^ 11 ^ 6 in the range 
M G [7.20 M Q , 7.60 M©]. We use flat priors on <j) c and 
t c over the range < (f>o < 2ir and —0.1 s < t c < 0.1 s, 
respectively. The prior on M. corresponds to the Jeffreys 
prior for the waveform family described by (2.5 1 [10 . 



For the mock data set we ran a MCMC comprising 
five parallel Markov chains in order to compute the PDF 
p(9 = (A4, r), D, t c , (/> c )\d,H). The limits of integration 
of the likelihood function, (2.2 1, are fixed to / m ; n = 40 Hz 
Hz, / max = 512 Hz. To extract the posterior samples 
from the raw MCMC output we discard the first 10, 000 
samples as burn-in. 

We measure the convergence of the parallel chains us- 
ing the Gelman- Rubin R-statistic [14]. For well con- 
verged chains this should be close to R = 1 and we regard 
the MCMC to be complete once R < 1.001 for all param- 
eters. 



SVD Setup 

The input to the SVD is a set of whitened time domain 
waveforms The frequency-domain SPA waveforms 
are whitened in the frequency domain with the PSD and 
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transformed into the time-domain for interpolation. By 
carrying out the interpolation in the time domain, we 
show that the technique can be applied to time-domain 
waveform families, though in this example we could work 
directly in the frequency-domain. Time-domain wave- 
forms are typically computationally expensive for param- 
eter estimation (see Sec.|V]), so this approach allows us 
to assess the computational savings associated with in- 
terpolating them. It is also consistent with the work in 
[U|9], where time-domain waveforms were interpolated. 

We ensure that all templates are of the same length, 
equal to the next highest power of two of the longest time- 
domain waveform in the set, which in our case is 2 s. For 
the proof of principle we apply the SVD to a small patch 
in Ai — r] space bounding the signal parameters. This 
region is set by the prior range on M. and r\ described 
above, chosen to be broad enough so that, for our data 
set, there is no likelihood support near the boundaries. 

The number of computations for the SVD of a N x L 
matrix with N < L scales like 0(LN 2 ). For the pur- 
poses of constructing the SVD we have found it effi- 
cient to split the M — rj space into four equally sized 
patches, with a separate SVD applied to each patch. 
The density of waveforms in each patches' bank is cho- 
sen such that the normalized inner-product between non- 
interpolated waveforms and interpolated waveforms gen- 
crated on an evenly spaced grid in each patch is at least 
99.9%. Such normalized mismatches of < 0.001 between 
interpolated and non-interpolated waveforms should en- 
sure that parameter-estimation accuracy is not compro- 
mised as long as the the signal-to-noise ratio does not ex- 
ceed ~ 20 (so that twice the mismatch times the square 
of the SNR is less than unity [15l US]), although pa- 
rameter estimation could remain accurate at even higher 
SNRs. For the mass space in this example, we find that a 
(15+1) x (15+1) grid of template waveforms in each patch 
is sufficient for the required accuracy. We add one to the 
grid in each dimension because the interpolation can be 
unpredictable at the boundaries of the space. The patch- 
ing is shown in Fig. [I] Each waveform in the template 
bank is generated at a fiducial distance of 1 Mpc. The 
truncated-SVD of the template bank in each patch uses 
20 basis waveforms. 

One subtlety of the implementation is that the interpo- 
lated waveforms are whitened time-domain filters. The 



likelihood function (2.2 1 is computed using un- whitened 



frequency series which in our example correspond to the 
Fourier transform of a complex time series. To recover 
the appropriate frequency series we first have to de- 
whiten the interpolated waveform by Fourier transform- 
ing the interpolated waveform and multiplying by the 
inverse of the PSD. While this process of de-whitening 
does not affect the normalized inner-product between in- 
terpolated and non-interpolated waveforms evaluated at 
the same point in parameter space, we have made no at- 
tempt to ensure the accuracy of the overall amplitude. 
We note that this is likely to have an impact on distance 
estimates. 



We further note that we do not use normalized wave- 
forms as input to the SVD as was done in [H [9] . Over 
large areas of the mass-space it is necessary to normalize 
the waveforms because there may be significant differ- 
ences in the power of waveforms across the space and 
so the SVD can disproportionally weight waveforms with 
more power if they are not normalized. We are not af- 
fected by this issue because all the waveforms in our ex- 
ample have roughly the same power due to the limited 
extent of the template bank in mass space. In general one 
would use normalized waveforms for SVD over a larger 
parameter space; however, to recover the overall ampli- 
tude of the interpolated waveforms the normalization co- 
efficients would themselves have to be interpolated. 

Below we compare the results of parameter estimation 
using interpolated and non-interpolated waveforms. 
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FIG. 1. Patching of the parameter space. The space is 
patched into four regions. The SVD is applied to each patch 
(shown in red) which is padded at the boundary by one wave- 
form and so contains a (15+1) x (15+1) grid of template wave- 
forms. The SVD of the template banks consists of 20 basis 
templates in each patch. 



V. RESULTS: COMPARISON OF PARAMETER 
ESTIMATES USING INTERPOLATED AND 
NON-INTERPOLATED WAVEFORMS 

The marginalized PDFs for complete MCMC runs 
using non-interpolated and interpolated waveforms are 
shown in Fig. [2j We have omitted the marginalized one- 
dimensional PDF of the coalescence phase (f> c as it of lit- 
tle physical interest. Each run required around 1.5 x 10 6 
waveform computations. The mean posterior values of 
the distributions along with the signal values are shown 
in Table |U 

The chirp-mass distribution computed using interpo- 
lated waveforms is clearly biased. This is corroborated by 
a two-sample K-S test which reveals that the two sets of 
samples are not consistent with arising from the same dis- 
tribution with overwhelming odds. Nevertheless, the sys- 
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Distance (Mpc) time (s) 



FIG. 2. Marginalized PDFs, ( |2.4[ |, produced using non-interpolated waveforms (blue or dark grey) and interpolated waveforms 
(red or light grey). The signal value is shown as the dashed red vertical line. 



Par am 


Mean posterior 
value (interpo- 
lated SPA) 


Mean posterior 
value (SPA) 


Signal Value 


M (M ) 


7.472 (2.5xl(T 2 ) 


7.467 (2.5xl0 -2 ) 


7.450 


V 


0.2457 
(7.1xl0~ 3 ) 


0.2457 
(7.2xl0~ 3 ) 


0.2473 


D (Mpc) 


32.39 (2.11) 


32.40 (2.13) 


33.00 


to (s) 


1.0 x 10~ 3 
(4xHT 4 ) 


1.0 x 10~ 3 
(4xl0~ 4 ) 






TABLE I. Maximum likelihood parameter estimates (and 
standard deviations) of the marginalized PDFs using inter- 
polated and non- interpolated waveforms (Fig.[2|. 



could be increased by, for example, using a higher den- 
sity template bank or using normalized waveforms as in- 
put to the SVD. In general the required accuracy can be 
estimated from the detection trigger SNR [T5] . 

The two-sample K-S test marginally fails for the 
coalescence-time distribution, but there is no evidence 
of a systematic bias on the scale of statistical measure- 
ment errors. We find that the sets of posterior samples 
for the other two PDFs in Fig. [2] symmetric mass ra- 
tio and distance, are consistent with arising from same 
distribution as quantified by the K-S test. 



Computational cost of template waveform 
generation 



tematic bias in the chirp mass is a factor of four smaller 
than the statistical measurement uncertainty. Thus, we 
pass a commonly-used threshold for sufficient waveform- 
model accuracy [e.g., [TB]. We note that the accuracy 



Two commonly-used time-domain waveform families 
which are relevant for parameter estimation are the 
inspiral-only post-Newtonian approximant Taylor T4 [12] 
and the effective-one-body family calibrated to numerical 
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relativity (EOBNR, [T2]) that includes inspiral, merger, 
and ringdown phases. The latter are typically more com- 
putationally intensive. 

For technical reasons, our comparison uses interpo- 
lated waveforms which are based on the SPA approxima- 
tion rather than time-domain waveform families. This is 
inconsequential because the SVD procedure is the same 
regardless of the form of the input waveforms. We have 
observed that applying the truncated SVD to Taylor T4 
time-domain templates in the mass space used to con- 
struct the SVD in Sec. |IV| yields the same number of 
basis vectors as when applying it to inverse-FFT'd SPA 
templates. The computational cost of interpolation will 
therefore be identical. 

Our measure to compare the computational costs of 
interpolated, TaylorT4 and EOBNR waveforms is the 
time it takes to compute a single interpolated waveform. 
While this does not estimate the theoretical minimum 
number of FLOPs of the process, and is also hardware 
dependent, it does provide a useful heuristic for com- 
paring the relative speed of each waveform family. Re- 
call that the interpolated waveforms are a time-domain 
approximant and hence the comparison is to determine 
the computational savings for time-domain waveforms. 
We restrict our comparison to waveforms generated in 
the mass space in Fig. [I] The length of TaylorT4 and 
EOBNR waveforms will in general depend on the spe- 
cific source masses. For a fair comparison we compare 
waveforms which have approximately the same number 
of data points. Because EOBNR must be generated at a 
sampling rate of 4096 Hz, we ensure that the interpolated 
and TaylorT4 waveforms are sampled at this frequency. 
Each waveform is approximately 2 s in duration. 

The results of the comparison are shown in Table [TT| 
For reference we also show the computational time of 
standard SPA waveforms. We find that on average, the 
interpolated waveforms are ten times faster to generate 
than Taylor T4 and fifteen times faster than EOBNR, a 
significant increase in computational efficiency. However, 
for the waveform parameters considered here, inspiral- 
only waveforms could be generated at lower sampling 
rates than the 4096 Hz required for EOBNR waveforms; 
therefore, the cost of constructing interpolated or Tay- 
lorT4 waveforms can be around four times smaller rela- 
tive to EOBNR than the values quoted in Table [H] 

We also estimate the cost of pre-computing the SVD 
interpolation. We have previously noted that the com- 
putational cost of an SVD of an N x L matrix with 
(N < L) scales like 0(N 2 L). One also needs to com- 
pute the N x L matrix of template waveforms as input 
to the SVD. The cost of computing a waveform of length 
L is typically O(L), possibly with a very large pre- factor. 
Thus, the total cost of pre-computing interpolation co- 
efficients will be less than 0(N 2 ) times the cost of an 
individual waveform computation. For instance, in our 
example, N — 16 x 16 = 256, so interpolation can reduce 
overall MCMC costs for any time-domain waveform tem- 
plates by an order of magnitude or more when the typical 



Waveform Family 


Computational Time (T) 


SPA 


0.2 


Interpolated 


1 


TaylorT4 


10 


EOBNR 


15 



TABLE II. Computational time of template waveform gen- 
eration in units of computational time of interpolated wave- 
forms, T. EOBNR, Taylor T4 and interpolated waveform fam- 
ilies are generated at a sampling rate of 4096 Hz and have a 
duration of 2 s. The interpolated waveforms consist of 20 pre- 
computed basis vectors. SPA waveforms are generated in the 
frequency domain; to ensure the SPA waveforms contain the 
same number of sample points they are generated at a sam- 
pling frequency A/ = 1/2 Hz and have a maximum frequency 
of 2048 Hz. 



MCMC chain length of > 10 6 samples is taken into ac- 
count. 



VI. CONCLUSION AND DISCUSSION 

We have provided a proof of principle that interpo- 
lated waveforms can be used for parameter estimation 
without unacceptable loss of accuracy. Our example was 
restricted to a five-dimensional search over the source 
chirp mass M. and symmetric mass ratio ?y, the dis- 
tance to the source D and the time and phase at co- 
alescence t c and 4> c . We further restricted the prior 
ranges on M and r\ to M G [7.20 M Q , 7.60 M Q ] and 
7] G [0.175, 0.250], respectively. The systematic biases 
observed when using interpolated waveforms are demon- 
strated to be smaller than statistical measurement uncer- 
tainties. Thus, SVD-interpolated waveforms satisfy the 
stringent waveform-model accuracy criteria imposed by 
parameter-estimation requirements. 

The relative computational times of generating in- 
terpolated waveforms and time-domain TaylorT4 and 
EOBNR waveforms are shown in Table [H] Interpolated 
waveforms can be generated at around an order of mag- 
nitude more cheaply than TaylorT4 or EOBNR. This 
suggests that the computational cost of parameter esti- 
mation can be significantly reduced by employing SVD- 
interpolated waveforms for likelihood computation when 
the latter is dominated by the cost of waveform genera- 
tion. 

In order for interpolated templates to be viable for pa- 
rameter estimation pipelines we need to apply the SVD- 
interpolation technique to a significantly larger region of 
the CBC mass space than in the example considered here. 
Searches of gravitational waves from low-mass systems 
look for binaries with a maximum total mass of 25 M 
and a minimum component mass of 1 Mq [17] and high 
mass searches target binaries with total masses between 
25 Mq and 100 M [T5]. To be able to apply our parame- 
ter estimation technique to triggers from such searches in 
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a single step, without first determining the more limited 
mass region where there is significant likelihood support, 
we will need to efficiently patch the parameter space over 
a large mass range so that the computational cost of gen- 
crating the SVD can be minimized. A necessary condi- 
tion for setting up the SVD in all patches is that its com- 
putational cost, plus the cost of generating interpolated 
waveforms, must be less than the cost of performing the 
parameter estimation with non-interpolated waveforms. 
This will be the subject of future work. 

Furthermore, we have to be able to extend the SVD to 
generic waveform families. Particularly interesting is the 
potential to extend the technique to EOBNR waveforms, 
which are currently expensive to generate, and waveform 
families which describe binaries with spinning compo- 
nents. The latter class of waveforms have an intrinsic 
parameter space with up to six more independent param- 
eters (two spin vectors) and the current technique of in- 
terpolation within the intrinsic parameter space of wave- 
forms may become costly in large-dimensional spaces. 
However, it is interesting to consider the potential to ap- 
ply the technique to spin-aligned/anti-aligned waveforms 
[e.g., 19, 20J as this class of waveforms have only one extra 



parameter, the reduced spin of the binary. The analysis 
of data from advanced LIGO and Virgo detectors, which 
may have lower-frequency cutoffs close to 10 Hz [J] , will 
require template waveforms that are several minutes in 
duration. This technique is likely to be highly relevant 
to parameter estimation in that context. 



ACKNOWLEDGMENTS 

We would like to thank Tyson Littenberg for reading 
through a draft of the manuscript and providing us with 
useful comments. DK is supported from the Max Planck 
Gesellschaft. KC is supported by the National Science 
and Engineering Research Council of Canada. RJES ac- 
knowledges support from a Perimeter Institute visiting 
graduate fellowship. Research at Perimeter Institute is 
supported through Industry Canada and by the Province 
of Ontario through the Ministry of Research & Innova- 
tion. This document has LIGO document number ligo- 
|pl200l36| 



[1] T. G. F. Li, W. D. Pozzo, S. Vitale, C. V. D. Broeck, 

M. Agathos, J. Veitch, K. Grover, T. Sidery, R. Sturani, 

and A. Vecchio, Phys. Rev. D 85, 082003 (2012). 
[2] I. Mandel and R. O'Shaughnessy, Class. Quant. Grav. 

27, 114007 (2010). 
[3] F. Bulik and K. Belczyski, ApJ. 589, L37 (2003). 
[4] G. Harry and The LIGO Scientific Collaboration, Class. 

Quant, Grav. 27, 084006 (2010). 
[5] Fhe Virgo Collaboration, 

Advanced virgo baseline design note vir027a0 (2009), 

URL https : //tds . ego-gw . it/itf /tds/f ile .php? 

callFil e=VIR- 0027A- 09 . pdf] 
[6] J. Abadie and The LIGO Scientific Collaboration, Class. 

Quant. Grav. 27, 173001 (2010). 
[7] V. Raymond, M. V. van der Sluys, I. Mandel, 

V. Kalogera, C. Rover, and N. Christensen, Class. Quant. 

Grav. 26, 114007 (2009). 
[8] K. Cannon, A. Chapman, C. Hanna, D. Keppel, A. C. 

Searle, and A. J. Weinstein, Phys. Rev. D 82, 044025 

(2010). 

[9] K. Cannon, C. Hanna, and D. Keppel, Phys. Rev. D 85, 
081504 (2012). 



[10] J. Veitch and A. Vecchio, Phys Rev. D 81, 062003 (2010). 
[11] M. van der Sluys, I. Mandel, V. Raymond, V. Kalogera, 

C. Roever, and N. Christensen, Class. Quant. Grav. 26, 

204010 (2009). 

[12] A. Buonanno, B. R. Iyer, E. Ochsner, Y. Pan, and B. S. 

Sathyaprakash, Phys. Rev. D 80, 084043 (2009). 
[13] D. Sigg and the LIGO Scientific Collaboration, Class. 

Quant. Grav. 25, 114041 (2008). 
[14] A. Gelman and D. B. Rubin, Statistical Science 7, 4, 

457-472 (1992) 7, 457 (1992). 
[15] L. Lindblom, B. J. Owen, and D. A. Brown, Phys. Rev. 

D 78, 124020 (2008). 
[16] F. Omne, Class. Quant. Grav. 29, 124002 (2012). 
[17] J. Abadie et al. (LIGO Scientific Collaboration and Virgo 

Collaboration), Phys. Rev. D 85, 082002 (2012). 
[18] J. Aasi et al. (LIGO Scientific Collaboration and Virgo 

Collaboration), arXiv: 1209.6533 (2012). 
[19] P. Ajith, Phys. Rev. D 84, 084037 (2011). 
[20] A. Faracchini, Y. Pan, A. Buonanno, E. Barausse, 

M. Boyle, F. Chu, G. Lovelace, H. P. Pfeiffer, and M. A. 

Scheel, Phys. Rev. D 86, 024011 (2012), 1202.0790. 



